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Abstract 

A quantum wave packet treatment of neutrino and neutral K and 
B meson oscillations is presented which incorporates the recoil particle 
in the production process, and includes the effect of the localization 
and lifetime of the source assumed to be a resonance or unstable parti- 
cle. This approach removes the ambiguities in the conventional single 
particle treatment of these oscillations, and elucidates the role of quan- 
tum correlations with the recoil particle. A fundamental connection 
between the stochastic decay time of the source and the space-time 
coordinates of the correlated final state particles is derived. 

A proper description of neutral K and B meson oscillations [|TJ pj and 
neutrino oscillations [H-pl requires that the familiar superposition of states 



with definite mass be represented by coherent wave packets || - flO|| . How- 
ever, the conventional single particle treatment of these oscillations leads to 
ambiguities which have lead to debates whether the momentum || - M or the 



energy |TT| -fl4| remains unchanged for different mass eigenstates, although 



the resultant transition probabilities are the same in both cases. Further- 
more, such descriptions leave unanswered the fundamental question how the 
properties of a single particle wave packet are determined by the nature of 
the production process. In this note we resolve these problems by consider- 
ing a a wave packet which includes a single recoil particle produced by the 
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decay of a resonance or unstable particle, and incorporates both energy and 
momentum conservation [15| -|18|| instead of resorting to conventional ad-hoc 



assumptions. Including final state correlations raises some interesting new- 
questions because these correlations take place between space-like separated 
events, a subtle problem in quantum mechanics discussed a long time ago 
by Einstein, Podolsky and Rosen (EPR) [fL9[] , and considered in connection 
with the production of a neutral KK or BB pair in |20| - [22]. In this case 



oscillations can be observed as EPR correlations between neutral mesons of 
fixed flavor or their decay products. It has been claimed in |23| - |26| that 
oscillations of the recoil particle can also be observed even if this particle has 
a fixed mass. In the case of neutrinos produced in the decay 7r — > /i + v where 
the recoil is a charged lepton, this would greatly facilitate experiments. How- 
ever, Lowe et al. ||27||-||28|| have argued that only a certain traveling pattern 
of oscillations in the recoil particle coordinates is observable, while Dolgov 



et al. |18| concluded that such oscillations could be observed only as EPR 
correlations provided that the detection of the recoil lepton is related to a 
neutrino of fixed flavor. 

These problems are examined here in a quantum mechanical wave packet 
treatment of neutrino and neutral meson oscillations which includes the recoil 
particle in an entangled state. We make explicit certain assumptions con- 
cerning correlations between two particles which have space-like separations. 
An important feature in this approach is that a correlated wave packet can 
incorporate the effect of both localization and finite lifetime of the source, 
which is assumed here to be a resonance or unstable particle. It will be shown 
that these properties explain why the propagation of the particles is confined 
near classical trajectories. While it has been recognized that classical mo- 
tion must be combined with wave properties and interference effects for an 
understanding of the oscillation phenomena, in current discussions classical 



trajectories have not been introduced in a self-consistent manner [jig] , p3 | - 
[ PH . A novel property of our wave packet is that it can incorporate the space- 
time coordinates of the decay point of the initial unstable state which can 
be observed and provides a reference point for the oscillations. Our results 
averaged over unobserved recoil particle coordinates, Eqs. |28| and |33|, differ 



from those of Dolgov et. al. ||18|| , who assumed that the source and decay 
particles satisfy exact classical equations motion in violaton of the principles 
of quantum mechanics. 

For simplicity we consider the theory in one dimension where all the 
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processes are collinear, and assume there are only two mass eigenstates \a > 
and | b > which is adequate for our purposes. Then the transition amplitude 
to some final state \g > is 

A oc cos{6) < g\a > ip a + sin(6) < g\b > ipb (1) 

where 9 is the mixing angle for an initial state of definite flavor or strangeness, 
and ip a and ipb are wavefunctions associated with the different mass eigen- 
states. We obtain these wavefunctions by time-dependent perturbation the- 
ory, assuming that the initial state is a resonance or unstable state of mass 
M and width T. The wavefunction for this state for t > is 

Mx, t) = J dpf(p)exp[ip(x - x s ) - (iE p + MT/2E p )(t - t s )} (2) 

where t s is the time at which this state is created as a wave packet centered 
at x s , f(p) is the amplitude associated with a momentum distribution p in 
the initial state with corresponding energy E p = ^p 2 + M 2 . In the following 
discussion we set for convenience x s , t s at the origin of our space-time coordi- 
nate system , but it should be remembered that in practice these coordinates 
are not known precisely and will be included in some of our main results. We 
assume that this amplitude has a sharp maximum at p — p, and expanding 
E p to first order in p — p we obtain 

■ip (x,t) = exp[i(px — Ept)] exp(—Mrt/2Ep) g(x — vt) (3) 

where the envelope of the wave packet is given directly by the wavefunction 
of the source at t — 

g(x) — exp(— ipx)ip(x, 0) (4) 
In the Wigner-Weisskopf approximation we obtain for t > 



where N = (l/2ir)JTM\vi2\/ E p , t> 12 is the mean relative velocity, and E x = 



p 2 + m 2 and Ei = yp 2 . + m 2 , are the relativistic energies of the two cor- 
related particles in the final state with masses mi and m 2 which can have 
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different values for the eigenstates labeled a an b. Conservation of total 
momentum in the production process implies that 



P = Pl+P2 



(6) 



A shortcoming of this representation for ip is that the state of both par- 
ticles is given at the same time t, while in practice these particles can be 
detected at different times t% and t 2 . However, since these particles are not 
interacting, the subsequent time evolution of the wavefunction can be deter- 
mined by their respective free particle Hamiltonians Hi and H 2 . Hence 

ip(x 1 ,x 2 ,ti,t 2 ) = exp(-iHi(t! -t) -iH 2 {t 2 -t))if>(x ll X2,t) 1 (7) 

and the required wavefunction ip(xi,x 2 ,ti,t 2 ) (footnote 1) is obtained by 
replacing the factor (Ei + E 2 )t in Eq. |5] by E x t\ + E 2 t 2 . If these particles 
are unstable, as is the case for neutral kaons, an additional factor in the 
integrand of Eq. [| is required, of the form 



where T\ and r 2 are the particle lifetimes. The assumption that the initial 
state is a resonance or unstable particle of width T implies an uncertainty 
in the total energy which plays an essential role in understanding how the 
decay particles are confined to classical trajectories. Nevertheless, we can 
define mean momenta p\ and p 2 associated with the mean total momentum 
p by the requirement that for these special values of momenta, the energy 
conservation relation 



is satisfied exactly. 

We carry out the integrations in Eq. [| approximately, by expanding the 
momenta p\ and p 2 around these mean values p% and p 2 , obtained as solutions 
of the momentum- energy conservation equations, Eqs. || and to first order 
in p — p and E — Ep, where E = E\ + E 2 . Second order terms contribute to 
the dispersion of the wave packet which we neglect here. We have 




(8) 



Ep = Ei + E 2 . 



(9) 



Pi = Pi + 



(E-E p )-v 2 (p-p) 



(10) 
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and 

(E-Ep) -vi{p-p) 
V2=P2 = , (11) 

where V\ = px/Ei, v 2 = P2IE2 are the classical or group velocities, and 
V12 = Vi — v 2 is the relative velocity of particles 1 and 2. Likewise 

E 1 = E 1 +v 1 (pi-p 1 ), (12) 

and 

E 2 = E 2 + v 2 {p 2 - ih)- (13) 

Changing the coordinates in the integrand of Eq. |5| to the variables p and 
E, we obtain for t\ 2 > 

ip(x 1 ,x 2 , t u t 2 ) = N'exp[i<j> 12 ) exp(-MTt 12 /2E p )g(z 12 ), (14) 



where N' = —iJTM/EpVi2 is a constant, and 

0i2 =pixi +P2X2 -Eih - E 2 t 2 , (15) 

t 12 = {Ax 2 -Ax 1 )/v 12 , (16) 

z l2 = (E 1 Ax 1 + E 2 Ax 2 )/E p , (17) 

where Axi = Xi — v{ti for i = 1,2. For t± = t 2 = t we have t%2 = t — [x\ — 
x 2 )/v\ 2 and z\ 2 = x — vt, where x = (E 1 xi+E 2 x 2 )/Ep is the center of mass of 
the two final state particles. Hence Zi 2 is the deviation of the center of mass 
from classical motion, and t\ 2 = td corresponds to the stochastic time td at 
which the particle pair is created at x\ = x 2 = Xd, where 4>i2 — P%d — Ept d and 
z n — Xd — vtd- At such a point the final state wavefunction is proportional to 
the initial wavefunction at x — Xd, t = td, and it is independent of the mass 
of the decay particles. Consequently this wavefunction satisfies the initial 
condition that the flavor of the state be independent of Xd, td, an important 
result which could not be imposed ab initio. As expected, for t — £ 12 < 0, 
this wavefunction vanishes. 

In our approximation, the initial function, Eq. determines directly the 
envelope g(zi 2 ) of the wave packet of the two final state particles. Hence the 
probability for finding these particles at x%,ti and x 2 ,t 2 is given by 

^-g 2 (z 12 ) exp(-TMt 12 /E p )dx 1 dx2 (18) 

EpVi2 
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It can be readily verified that the probability for the creation of these two 
particles at xj, td is equal to the probability that the source decays at this 
same space-time point. If this decay is not measured or constrained by the 



environment, then Eq. [18] can be applied directly to calculate probabilities 
or averages over the recoil variable coordinates. However if is observed 
then for times ti > td and t 2 > td the probability for finding the particles at 
Xi,x 2 is obtained by setting t 12 = td in Eq. [18] (footnote 2), and x\ — x 2 = 
vi(h - t d ) - v 2 (h - t d ), or 

A d x 1 = A d x 2 = z 12 (19) 

where AdXi = Xi — vtd — v%{ti — td) is the deviation from classical motion of 
the decay particles. In this case, the measure zi 2 for this deviation satisfies 
he same distribution as the deviation from classical motion of the initial 
state, and therefore our analysis shows that particles associated with the 
decay process are confined to move along classical trajectories with the same 
degree of localization as the source. 

We now assume that the two states labeled a and b correspond to parti- 
cles 1 and 2 with small mass differences, 5m^ = mf a — mf b1 and calculate the 
corresponding differences in the mean momentum pi and p 2 and correspond- 
ing energies E\ and E 2 from the energy-momentum conservation laws, Eqs. 
H and H We have 

Spi + 5p 2 = 0, (20) 

and 

6E t + 5E 2 = 0, (21) 

where to first order in 5m? 



5E t = + (22) 



2Ei 



Solving these equations we obtain 



5 Pl = -6p 2 = -±( 5 4+ 5 4). (23) 
y y v l2 y 2E l 2E 2 ! y 1 

where V\ 2 = V\ — v 2 is the relative velocity. These relations differ from the 
result obtained with conventional kinematics assumptions that different mass 



states have either the same momentum p] - H or the same energy [11] -[14 
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The oscillation term which concerns us here appears in the calculation of 
the interference term in the transition probability \A\ 2 , where A is given by 
Eq. [TJ, and is proportional to 

Reatyl^b = cos(4>)g(z l2A )g{zi2,b)exp{-MT{ti 2) a + t 12 , b )/2E p ) (24) 

where we have ignored factors which depend on the lifetime of the final state 
particles, Eq. ||. The phase difference = 4>u,a — 0i2,f> = fipiXi + 8p 2 x 2 — 
5Eiti — 5E 2 t 2 is invariant under Lorentz transformations, and according to 



the energy-momentum conservation laws, Eqs. |2(J and |2l], it can be written 
in the form 

<j) = 5p 1 (x 1 -x 2 )-5E 1 (t 1 -t 2 ). (25) 

This form shows that the phase difference depends only on the relative time 
coordinate of the final state particles and, therefore, that it is independent 
of the initial decay point td. However the role of this decay time appears 
when we consider the effect of the wave packet envelope in the case that this 
coordinate is measured, constraining each of the decay particles to move near 
its classical trajectories. For this purpose we apply Eq. |22| to write in the 
equivalent form 

= 8p 1 (Ax 1 - Ax 2 ) - -^±ti - 2^t 2 , (26) 



This expression for is similar to results given in JL8] and p8[ , but our 
correlated wave packet now allows us to interpret and evaluate properly the 
contribution of the first term in which does not appear in the convential 



single particle formulation for this phase. Substituting Eq. |23] for 5p\ and 
substituting Ax 2 — Ax± = vi 2 ti 2 , we obtain finally 

- Sl %,-^-^-t 12 ), (27) 



2E 1 y 1 2E 2 

If the recoil particle has a fixed mass, i.e. 8m\ = 0, this form for is 
equal to the conventional single particle result with mass eigenstates of the 
same momentum provided that we identify t\ 2 with the decay time td of the 
initial state as we have done previously, see footnote 2. In this case this 
phase is independent of the recoil particle coordinates contrary to assertions 
^]-|f26|. However, if the decay time td is not measured directly, then in 



m 
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principle it could be determined from a coincidence measurement of the recoil 
coordinate x^- If we set t\ = £2, and substitute t\ — tyi = {x\ — x-i)jv\i in 
Eq. ^7], we obtain EPR-like oscillations in the relative coordinates of the two 
final state particles which can only be observed if the flavor is also determined 
[ |l8|j . In the case of neutral meson pairs in the final state, e. g. KK produced 
in <fr decay, b~m\ = —8m\, and the time coordinates of both particles appear 



in Eq. 27 even if the decay time td has been determined. A similar result 



was obtained in |^TJ and |22] by assuming that Spi = Sp 2 = 0, although this 



kinematical condition is not justified. 

If neither the decay time td nor the coordinates X2,t 2 of the recoil par- 
ticle are observed, we must integrate the interference term of the transition 
probability, Eq. ^ over the decay time td, or what amounts to an equivalent 
procedure, over the position coordinate x 2 of the unobserved recoil particle. 
Assuming that there are no enviromental constraints on the possible range 
of these variable, we apply eq. [18] to obtain the average of cos(0). Neglecting 
the mass difference in the envelopes of the wave packet, which leads to a 
finite coherence length, we obtain 

< cos(<P) >= R cos (~(ti - t s ) + -J^(t 2 - t s ) - 5), (28) 

where 

R= . 1 „ , (29) 
tan{ S) = f , (30) 

and 

This average is independent of the shape of the initial wave packet. The 
parameter £ gives a measure for the deviation from the conventional form, 
which corresponds to R = 1 and 5 = in Eq. |28|. For example, in particle 
reactions producing neutral mesons, T is of order several MeV, and for the 
B meson 5m is 3.1 x 10~ 4 eV , and about 100 times smaller for the K meson. 
Hence £ is of order 10~ 10 — 10 -12 , and the contribution from the first term 
which appears in Eq. |26| is essentially unobservable, contrary to expectations 
in ||28|| . For neutrinos produced in pion decay, with 5m 2 ~ lO -3 ^ 2 as found 
in |J, £ fa 10~ 3 can also be neglected in the analysis of the data. However, 
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this is not the case for 5m 2 of order 1 — 10 eV 2 as reported in M. Furthermore, 
the magnitude of £ is greater for the case that the neutrinos are produced in 
muon decay [|3| because the muon is lighter than the pion and has a longer 
life time, but this is a three body decay and our analysis is only approximatly 
valid. 

In practice it must be remembered that oscillations are observed in po- 



sition rather than in time coordinates. Setting ^ = (x{ — Axi)/v{ in Eq. |27 
for the phase we obtain an equivalent form 

,8m 2 5ml. 5m 2 . . 5m 2 ,. . . . 

* = ( %f + ^ )212 " ift (x > ~ mv ' ] ~ ^ ~ mn) - (32) 

This differs from the conventional form for 0, but if the decay coordinates 
are determined then tn = tj, Zxi = A^xi = A^X2, Eq. [19], and we recover 
the conventional form. Notice that in the rest frame of the source v = 0, 
and the dependence of the phase on t 12 vanishes because the source is not 
moving. Assuming a Gaussian distribution for the initial wave packet with 
a width a x , and averaging cos(4>) over Z\i and t\i we now obtain 

< cos(4>) >= R cos ( — —(xi — x s ) H —{xi — x s ) — 5), (33) 



2p! y ' 2p 



2 



R = ^tfi, (34) 



and 



where 



tan( 5) = £, (35) 
,5m 2 5ml. 

As expected from simple physical arguments, interference effects can occur 
in position measurements provided that the width a x of the wave packet is 
small compared to the oscillation length oc pi /5m 2 , or correspondingly that 
r] oc 5pi<j x « 1. The magnitude of a x is of the order of magnitude of the 
localization of a nuclear target, and it is further contracted by the Lorentz 
transformation due to the motion of the unstable initial particle so that in 
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practice 77 << 1. Moreover, in the rest frame of the source p = and 
consequently in this case £ = 0. 

In conclusion, we have shown that the transition probability for neutrino 
and neutral meson oscillations can be obtained from first principles by solv- 
ing the time dependent Schrodinger equation for the decay of an unstable 
source into a coherent superposition of correlated two particle eigenstates 
with different masses. We have obtain our results in a relativistically co- 
variant manner by applying well defined approximations without recourse to 
conventional ad-hoc assumptions which violate principles of quantum me- 
chanics, and have led to much confusion in the literature. We have shown 
that the width or lifetime of the source plays a crucial role in understanding 
this problem, and that quantum correlations between the final state particles 
relate the decay time of the source to the space-time coordinates of the these 
particles. Due to the mass difference the wave packets for different mass 
eigenstates have different group velocities and separate leading to a finite 
coherence length [fL5||, but this effect was neglected here. 
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Footnotes 



The justification for Eq. [7] is that after one of the two particles has 
been detected, i.e. t — t\ or t — t 2 , its state does not continue to evolve 
in time. The generalized wavefunction ip(xi, x 2 , ti,t 2 ), Eq. [TJ| can then 
be interpreted as the probability amplitude for correlated events which 
occur at the two different space-time points (xi,ti) and (x 2 ,t 2 ), and is 
equivalent to the amplitude method in ||21|| , and the formalism in [|l^ 



It can be shown that this procedure is equivalent to the "collapse" of 
the wavefunction language, which is the conventional description when 
measurements take place at different times, although it is preferable 
not to invoke this awkward language. 

2. For ti,t 2 > td-, the condition t 12 = {v\ti - v 2 t 2 — xi + x 2 )/v 12 = t d 
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connects the decay time td of the source to the space-time coordinates 
of the decay particles. It corresponds to the classical relation for the 
relative coordinates of these particles which can be understood on the 
physical grounds that these particles are created in a region of negligi- 
ble small spatial dimension, without violating the uncertainty principle 
because td is a stochastic variable. This relation has also been obtained 
in [18|], but under the invalid assumption that the source and decay par- 
ticles follow classical trajectories exactly. In the case that the decay 
can occur over a range of values < td < t max , one must take an aver- 
age over the probability distribution integrated over this range provide 
there are no measurements on the recoil particle which constrain the 
possible values of tu. In practice t max = d/v, where d is the distance 
between the target where the unstable particle is created, and a beam 
stop where nuclear reactions annihilate it. 
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